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>- ■ Two-nucleon scattering at intermediate energies of a few hundred MeV requires quite 
a few angular momentum states in order to achieve convergence of e.g. scattering ob- 
servables. This is even more true for the scattering of three or more nucleons upon each 
other. An alternative approach to the conventional one, which is based on angular mo- 
mentum decomposition, is to work directly with momentum vectors, specifically with the 
magnitudes of momenta and the angles between them. We formulate and numerically 
illustrate this alternative approach for the case of NN scattering using two realistic 
interaction models, the Argonne AV18 and the Bonn-B || potentials. The momentum 
vectors enter directly into the scattering equation, and the total spin of the two nucleons 
is treated in a helicity representation with respect to the relative momenta q of the two 
nucleons. 

v~j ■ The momentum-helicity states are given as 
rN ' , , 

\q;qSA) = |q> \qSA) = |q> R(q) £ C(--S ] m l m 2 k) 
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where R(q) = exp(—iS z (j)) exp(—iS y 9) is the rotation operator, S z , S y the components of 
S = |(<Ti + cr 2 ), and A the eigenvalue of the helicity operator S • q. Introducing parity 
and two-body isospin states \tm t ), the antisymmetrized two-nucleon state is given by 

|q; qSA; ty a = -L(i - v A-) s+t ) \t) |q; qSA)„ , (2) 

with the parity eigenvalues 77^ = ±1 and |q; qSA) n = ^(|q) + ^ |— q)) I^A). 

With these basis states we formulate the Lippmann-Schwinger (LS) equations for NN 
scattering. In the singlet case, S = 0, there is one single equation for each parity, 

2#«(q', q) = ^(q\ q) + \ J d 3 q" , 4')G (q")T f *(q", q) (3) 

Using rotational and parity invariance one finds in the triplet case, S = 1, a set of two 
coupled LS equations for each parity and each initial helicity state 
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:q, q) = vfM, q ) + 1 -J d 3 q" vz s M, ^)G ( q ")T^', q ) 

+\ J dY VfM, <f)Ga(q")T T (q", q). (4) 
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Because of the symmetry properties only two (A = 1, 0) of the original three initial helicity 
states need to be considered. Both LS equations are two-dimensional integral equations 
in two variables for the half-shell t-matrix and in three variables for the fully-off- shell 
t-matrix, namely two magnitudes of momenta and one angle. For explicit calculations 
we choose q = z, which allows together with the properties of the potential that the 
azimuthal dependence of q' and q can be factored out 

T A ^(q', q) = e iA ^Tfl(q', q, 6'). (5) 

As it is well known, rotational-, parity-, and time reversal invariance restricts any NN 
potential V to be formed out of six independent terms . We introduce the following set 
of linear independent operators, which are diagonal in our basis 



Q 1 = 1 
Q3 — S • q' S • q' 
Q 5 = (S • q'f (S • qf 



n 2 = s 2 

Q4 = S • q S • q 

Qq = S • q S • q (6) 



and represent the most general potential in a nonrelativistic Schrodinger equation as linear 
combination of them, 

(q'|V|q) = V(q',q) = £ Vi (q', q, 7 )0,. (7) 

i=l 

Here Vi(q', q, 7) are scalar functions of the magnitudes of the vectors q', q, and the angle 
7 = q' ■ q between them. It should be noted that for q pointing in z-direction one can 
use (q'SA' \zSA) = exp(iA(0' — 0))<if A ,(6*'), so that the azimuthal dependence factors out 
of V and thus also out of the T-matrix. The operators fl{ can be transformed to the 
standard Wolfenstein form [Bl with a linear transformation as shown in llifl. 
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Figure 2. | (++|T| )| 2 and | (++|T|++) | 2 in units of 10" 14 MeV" 4 

as a function of x = cos 6* for g =375 MeV/c. 

To demonstrate the feasibility of our formulation when applied to NN scattering we 
present numerical results for two different NN potentials of quite different character, the 
Bonn-B || and the AV18 @ potential models. Similarities and differences in the half-shell 
T-matrices TJ[^(q, q , 6) are displayed for T^ 1 (q, q , 6) for q = 375 MeV jc in Fig. 1. The 
amplitudes are identical on-shell and show deviations off-shell. The choice presented in 
Fig. 1 is quite representative for the off-shell differences between the two models. 

In order to connect with standard representations of the NN t-matrix as well as to 
calculate observables, we express T in terms of antisymmetric states 

|rir 2 mim 2 q) a = -^=(1 - P 12 ) |rir 2 mim 2 q) , (8) 

where Ti and m« are the magnetic isospin and spin quantum numbers and Pi 2 is the 
permutation operator of the two nucleons. In this representation, which we denote as 
'physical', we obtain for the on-shell T-matrix elements ( |q'| = |q|) 

.(wiXffS'l T |r 1 r 2 m 1 m 2 q) a = l e -^-^ £ C(~t; nr 2 ) 2 (l - ^(-) 5+t ) 

Sirt 

C{\\S- m'^Kp^S; m im2 A ) £ df^(ff)T^(q, q, 6'). (9) 

The quantity Tj[F£ (q, q, 6') can be related through straightforward but lengthy algebra to 
the standard partial wave T-matrix TJ,/* (q) as is shown in |I[ . Here we want to consider the 
on-the-energy-shell amplitudes \ a (riT2m' 1 m' 2 q q\T\TiT2mim 2 qo)a\ 2 = | (m' 1 m / 2 |T|mim 2 ) | 2 , 
m = ±|. If one takes rotational symmetry and parity invariance into account, one ends 
up with six independent amplitudes out of the 16 possible combinations. Two of them, 

|(+ + \T\ )| 2 and |(+ + \T\ + +)| 2 , are displayed in Fig. 2 as a function of cos# where 

6 is cm. angle for a projectile energy of 300 MeV (go=375 MeV/c). In addition to the 
result from the 3D calculation we show the partial wave sums for increasing j ma x- One 
can clearly see that at £"^6=300 MeV partial WclVGS lip to dbt ].6clSt Jmax — 12 are needed to 
obtain a converged result for every matrix element. 
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From the 'physical' T-matrix amplitudes we construct the Wolfenstein amplitudes |J. 
Once those are obtained it is straightforward to calculated the NN scattering observables 
|J. In order to unambiguously define the normalization we also give the expression for 
the spin-averaged differential cross section for nucleon species tit 2 as 

^ = ( 27r ) 4 (tt) \ E \a{riT 2 m' l m' 2 qq l \T\T l T2m l m 2 ^i a \ 2 (10) 

In Fig. 3 the differential cross section is given for .£^=300 MeV together with partial 
wave sums for increasing j m ax- It should be noted that the convergence of the partial 
wave sums is especially slow for the backward angle differential cross section. A sum up 
to j m ax=^G is needed to obtain convergence within 1 %. In Fig. 4 A y is displayed, which 
needs a partial wave sum up to j ma x=10 to reach a similar convergence. 
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Figure 3. % as a function of 9 for 
g =375 MeV/c. 



Figure 4. A y as a function of 9 for 
g =375 MeV/c. 



In summary, we formulated and numerically illustrated an alternative approach to NN 
scattering which works directly with momentum vectors. The spin of the two nucleons 
is treated in a helicity representation with respect to the relative momentum of the two 
nucleons. We would like to emphasize that the here developed scheme is algebraically 
quite simple to handle provided potentials are given in an operator form. This is e.g. 
the case for all interactions developed within a field theoretic frame work. This work is 
intended to serve as starting point towards treating three-nucleon scattering in a similar 
fashion. 
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